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1 Introduction 

The solution of the biobjective IRP is rather challenging, even for metaheuristics. We are still lacking 
a profound understanding of appropriate solution representations and effective neighborhood structures. 
Clearly, both the delivery volumes and the routing aspects of the alternatives need to be reflected in an 
encoding, and must be modified when searching by means of local search. Our work contributes to the 
better understanding of such solution representations. On the basis of an experimental investigation, the 
advantages and drawbacks of two encodings are studied and compared. 

2 Representations 

Due to the high complexity of the problem, many solution approaches introduce some simplifying as- 
sumptions. Some authors restrict quantities to deliver |2j or allow backlogging [IJ. Cycling visiting was 
also studied [S] and other authors limit routing possibilities [5]. 

To reduce the combinatorial complexity, we choose a specific replenishment strategy such that shipped 
amounts are not part of the searched variable set anymore; moreover, visits are periodic [4J. In our first 
representation, only the visiting frequency fi of each customer i has to be decided; shipping amounts are 
deduced from our replenishment strategy. Iir our new representation, the starting date Si of the periodic 
visiting sequence is adjustable; one delivery before the starting date is allowed if required to cover first 
consumptions of some customers. 

Our search algorithms will have to explore values of n or 2n variables, instead of up to Hn variables 
for other approaches {H is the size of the time horizon) . Such simplifications make us miss solutions and 
in particular some non-dominated ones. Nevertheless, we expect that the loss of freedom in quantities 
and visiting dates is balanced by the higher local search speed. 

3 Local Search and Evaluation 

Not surprisingly, solution approaches often employ heuristics and meta- heuristics [7]. Following in these 
footsteps, we began to experiment our representation with an heuristic. For the moment, our neighbor- 
hood operator is simple and we investigate representations and solution selection strategies. 

The first selection strategy uses reference points [S]. At each iteration, among all non dominated 
solutions, we select the R closest solutions to R evenly spread reference poiirts. Our second strategy is 
inspired from NSGA-II [3]. At each iteration, among all non dominated solutions, a subset of R solutions 
is randomly chosen. The probability that a solution is taken depeirds on its crowding distance. 

Neighbors are generated for each of the R alternatives then inserted among all other solutions and 
dominated solutions are removed. This process is repeated until the total number of evaluations reaches 
a limit chosen by the user. 

The local search adjusts starting dates and frequencies of the visits. From our short representations, 
the role of the evaluation phase is to compute holding costs and routing costs of each solution. Whereas 
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inventory costs are fast to compute exactly with our replenishment strategy, the routing cost evaluation 
is a NP-Hard problem and we must restrict ourselves to a rough value. 

The sum of inventory levels L* of all customers i over all dates t represents the overall holding cost: 

zi — XltLi 127=1 ^i- Evaluation of inventories takes Q{Hn). 

Concerning routing costs, we get them from the sum of the costs of the VRP solutions for all dates: 
Z2 = X^^i VRP((7*), where is the set of quantities to deliver to customers visited at the date t. 
The route planning is preliminary solved by the Clarke & Wright's heuristic, which allows for a fast 
construction of tours in 0(n^ logn). Then, a 2-Opt operator optimizes the tours independently in O(n^) 
in average. 



4 Experiments and Results 

We denote by ev the number of solutions evaluated since the beginning of the search. It is interesting 
to note that the number of modifications done on any initial solution during search is at most and 
that the average number of modifications applied per customer can not exceed 4^2^. In our previous 
approach, a customer was served for the first time just before going to stockout [4J. As described above, 
starting dates are now movable. We noticed that this new freedom for delivery dates is always beneficial 
when 4^td is constant with regard to n. This is a typical result on the left: 
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On the right, we notice that our random selection method reduces the gaps, but sometimes at the expense 
of the solution quality. 

In addition to those experiments, we defined the similarity between two solutions by considering either 
their values in the objective space or in the decision space. The purpose is to remove identical solutions. 
Better results are obtained if only the objective values are compared. 
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